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1 Introduction 

The well-known Hadamard method of descent ]2|, ^] consists of deriving solutions 
of wave equations in n variables from (known) solutions in more than n variables. 
To quote Hadamard, (the method) "consists in noticing that he who can do 
more can do less" . It is the purpose of this paper to exhibit cases in which he 
who can do less can do more — that sometimes one can build the propagator 
cos [\/A 2 + B 2 t) out of cos(Ai) and cos(Bt). 

In Section 2 we consider the commutative case and represent 



cos ( yj A\ + ■ ■ ■ + A 2 t tj in terms of an integral involving cos(Aiti) (Theorem 2.1). 
We distinguish between the cases n even and n odd. While it is true that once 
you see the formula ((2.1) or (2.2)), you may verify it directly (and this is the 
method used for the case of n even), we prefer to use for the case n odd a 
method which shows how the trigonometric-hyperbolic case is derived from the 
elementary exponential (parabolic) case exp (— p(A\ + • • • + A%)) through "in- 
version" of a transmutation formula . The case n even could also be reduced 
to the case n odd via a descent argument. Note that (2.1) and (2.2) appear to 
be new even for the case of numbers (scalar operators). The cases n = 2, 3 are 
rather simple: 




cos 



A 2 + B 2 t 



) 




cos(tcuiA) cos(tw2-B) 



y/l - (iU 2 + LU 2 ) 



du>idu>2 



(1.1) 





(1.2) 



t / cos(tcoiA) cos(iw2-B) cos(tLU3C)duj 
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As a simple illustration, we derive the well-known formulas for the solution 
of the initial value problem for the wave equation in R n from the elementary 
solution for the one-dimensional case (Example 4.1). Note that we do not use 
explicitly the spherical symmetry of A. We also solve the initial value problem 
for the Klein-Gordon equation. 

In Section 3 we derive formulas for cos (V A 2 + B 2 t) when A and B do not 
commute. That the expression should involve a limit follows from the Trotter 
product formula for e~^ A +B '. Formal "inversion" of the Trotter formula leads 
to a limit formula for cos (V A 2 + B 2 <) . This formula may be justified weakly 
(Proposition 3.1) or strongly on appropriate analytic vectors || (Theorem 3.2). 

One could express cos {^J A\ + ■ ■ ■ + A 2 tj for general n (Remark 3.4) but the 

resulting formulas are quite heavy. The formulas (such as (3.4)) are reminiscent 
of those appearing in path integrals (such as the Feynman-Kac formulas for 
parabolic or Schrodinger equations) . This similarity consists of the integrations 
being extended over balls of increasing dimensions; due to the hyperbolic nature 
of the problem, we have here in addition also differentiations to high orders. 

To illustrate the formula for the non-commutative case we compute the prop- 
agators of the harmonic oscillator 

d 2 u d 2 u 2 
dt 2 dx 2 

(Example 4.2) and of simple hypoelliptic sum of squares operators, such as 

d 2 u d 2 u 2 d 2u 
W = dx\ +Xl dxl 

and 

d 2 u . 

where A# is the Laplacian on the Heisenberg group (Example 4.3). The ex- 
pressions obtained are not very explicit. 

sin(V A' 2 + B' 2 t) 

Mote that the operator — \/a 2 +b° ' a more generally, the operator 

sin[y/A\ + -+Al t) 



y /A 2 1 + -+A$ 



may be represented by similar formulas, the only difference be- 



ing dropping the left-most from the corresponding formulas for 

cos (v 'A 2 + B 2 t) ^cos (t^ A\ + ■ ■ ■ + A 2 ^ ^ . The proofs are similar to those 

of (2.2) or of Theorem 3.1. 

I am very much indebted to V. Katsnelson for simulating discussions. 

2 The commutative case 

Theorem 2.1 Let H be a Hilbert space, and let the {possibly unbounded) self- 
adjoint operators Ax, ... , A n commute. 
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(i) // n is odd, n = 2m + 1, then 



coBltJAl + --> + A* = 



1 9/13 



2(27r) m dt\tdt 
x / cos{tw\Ai) ■ ■ ■ (tuo n A n )duj 



t2m-l 



(2.1) 



(here du> denotes the surface measure on S ) 
(ii) If n is even, n — 2m, then 



t^Al + .-. + Al 



1 ra 9/19 



,2m-l 



(2tt) dt\tdt 

COs(fcJlAl) • • • COs(tui2mA2m) 



dwi • • • du)2r, 



(2.2) 



Proof. Set L = VEILi A ? • Note that L 

is a non-negative self-adjoint oper- 
ator in H. To prove (i), note that by commutativity we have for p > that 



n« 

i=l 



The well-known transmutation formula (e.g. [Q) 



/4irp 



valid for B self-adjoint, implies that 



-B 2 P _ 



1 



npJo 



e *p cos(Bt)dt , 



e 4 p cos(Bt)dt , 



(2.3) 



(2.4) 



(2.4') 



Applying (2.4) to each Ai, we find that 



n« 



-A* iP _ 



e 4 p cos(j4iti)dti / e 4 p 008(^2*2)^2 



(47rp)"/ 2 

/'CO 

x / cos(A n t n )dt n 

e 4 p cos(A±tx) ■ ■ ■ cos(A n t n )dt 



1 



(47rp)™/ 2 
1 

(^P) n/2 J 



t n e 4 p / cos(tojiAi) ■ ■ ■ cos(tuj n A n )dujdt 



(2.5) 
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Setting B = L in (2.4') we see from (2.3) and (2.5) that 

1 f°° t 2 1 f°° i 2 

— e~-p cos(Lt)dt =- — -7TT / t n ~ x e~~^ 
^PJo (4^)™/2 J Q 



/ cos(tujiAi) ■ ■ ■ cos(tuj n A n 
Js™- 1 



)dudt 



so that 



where 



J e"*? cos(Lt)dt = 2{A7rp)m J e-^t 2m F(t; A u ..., A n )dt (2.6) 

F(t; A\, . . . , A n ) = I cos(tujiAi) cos(iw2^b) • • • cos(tLu n A n )cLo . 
Js 2 ™ 

We apply successive integrations by parts to eliminate the factor p~ m that 
appears in the right hand side of (2.6). Thus, 



1 



e--?t 2m F{t;A l ,...,A n )dt 
1 1 f°° d 



2?r (4irp) m - 1 J dt 
1 1 



e~i-At 2m - 1 F(f,A 1 ,...,A n )dt 



r°° t 2 r) 

i y e-^^- t [t 2m - 1 F(t;A 1 ,...,A n )]dt 



2ir (AttpY 

(t 2m - 1 F(t;A 1 ,...,A n )]dt 



i i r 

(2tt) 2 (4np)™-2j 
(27r) m i 



e 4 p 



di\tdt 



' d (\ d 

e 4p — 

dt\tdt 



m — 1 



[^-^(Mi,---,^)] <ft 



Substituting in (2.6) and applying the uniqueness theorem for the Laplace trans- 
form, we obtain (2.1). 

One could derive (2.2) from (2.1) by setting A 2m +i = (essentially a method 
of descent). We prefer an alternative method. Let ai, . . . , a 2m be scalars. Then 
(using multi-index notation) 



/ 



cos(wiaii) • • ■ cos(aj2 m a2mt) 



v/1 



du>i ■ ■ ■ du>- 



2m 



f2|a| a 2a f 



(2.7) 



(2a)! i|| w ||<i V 1 - II w || 2 



dwi • • • duj2 
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The integral on the right hand side of (2.7) is essentially the Dirichlet integral 
[|. In fact, 

,,,2a r w 2a 

duii ■ ■ ■ du>2 m — 2 2m / — =du>i ■ ■ • duj n . 



\\ 2 <1 || U} || 2 ' Jlu>0, \\lu\\ 2 <1 V 1- II W I 

Set Zi = w 2 , 1 < i < 2m. Then the integral on the right hand side of (2.7) is 
equal to 

Z T 2 "' z 2m m 2 (1 ~ Z l _ • ■ ' ~ Z 2m y 1,2 dzi ■ ■ ■dz 2m ■ 

z>0, «3m<l 

According to formulas (7.7.4) and (7.7.5) in ||, the value of this integral is 

r(ai + |)---r(tt 2m + i)r(i) 

r(|a|+m + i) 

Substituting in (2.7) and recalling the duplication formula for the Gamma func- 
tion _ 

1\ y^T{2k) 



1 ' k ' 2)~ 22fe-ir(fc) 



we see that 

COs(wiait) • • -COs(a; 2m a 2m i) 

|| 2 <1 v 71 " IMI 2 



du>i ■ ■ ■ du>[ 



2 m 



E 

a 

E 



(2a)! r(|o|+m + i) 

(-1)H ^l"la 2a 7r m r(2ax) ■ ■ ■ r(2a 2m )2 2 l a l+ 2m - 1 r(|a| + m) 
(2a)! 2 2 H" 2m r(ai) • • • r(a 2m )r(2|a| + 2m) 

(_l)l"lt2|a| a 2a 7r m 2 4m 1 

(2ai) • ■ • (2o 2m )r(ai) • • • I> 2m ) 2M+™(2m + 2\a\ - 1) 
1 



E 

a 

E 



(2m + 2|a| - 3) ■ • • 1 

(_ 1 )|a| t 2H a 2a 7r m 2 4m j j 

al 2 2m 2l Q l+ m (2m + 2|a| - l)(2m + 2|aj - 3) • ■ ■ 1 

(-l)l Q li 2 l a la 2a (27r) m 1 

2l Q la! (2m + 2|a] - l)(2m + 2|a| - 3) • • • 1 ' 
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Hence 



'/ 

J II, 



.,,„ | / cos(wiaii) • • • cos(u> 2m a2 m t) 

t I (LU)\ ■ ■ ■ (Md2m 

W |p<i ^1- |MI 2 

(-l)l"l(27r) m t 2|a|+2m-l a 2a 



(2.8) 



E 



2Ha! (2|a| + 2m - l)(2|a| + 2m - 3) • • • 1 ' 



But 



1 d 



t dt 



~ W2|a|+2 m -l = (2 | a | + 2m _ 1)t 2\ a \+2rn-^ ...It") p\a\+2 m -l 



i a 



1 9* 



m— 1 



= (2|a|+2m-l)---(2|a|+3)f 2|Q|+1 . 
Applying to (2.8), we get 

cos(wiait) • • • cos(co2 m a2mt) 



ld_ 

Idt 



m— 1 



= E 



'ii^ii 2 <i 

(-l)l"l(27r) m a 2a 



VHFF 

f2|a|+l 



(iWi • • • d0J2m 



2Ma! (2|a| + l)(2|a| - !)•••! 



and 



1 8 fid 



(27r) m dt\tdt 



m—l 



COs(^iflif) ■ ■ ■ COs(u) 2m a2mt) 
||u-P<l y/l- II W IP 



dwi • • • dw 2m 



L I Ola 



(_l)l«l a 2a^2|a| 

2l Q la!(2|a|-l)---l ' 
On the other hand, for every non-negative integer fc, 



(2.9) 



Inserting in (2.9) we conclude that 



1 d (Id 



(27r) m dt\tdt 



m—l 



i2m-l 



cos(wiait) • • • cos(cj 2m a 2m t) 



ll 2 <i 



dwi • • • dui 2m 



~ (-l)*t2*(o? + "- + aL) fc ^ (-l) fc (|«|t) 2fc 



2 fc fc! (2fe-l). 



E 

fe=i 



(2fc)! 



cos(|a|t) 



proving (2.2) for the case Ai = dj, i = 1, . . . , 2m. In general, let dE^(X) denote 
the spectral measure of Ai, 1 < i < 2m. The commutativity assumption implies 



G 



that 



cos tJAf 



Application of (2.2) (case of real numbers) to (2.10) implies that 



cos tJAf 



1 9/18 



(27r) m dt\tdt 

COs(wiAlt) • • • COs(w2mA2mi) 



2m — 1 



,2m-l 



„ ii 2 <i V 1 - II u I 

xdE^(X 1 )---dE^ m \X 2m ) 



dw\ ■ ■ ■ dw2m 



1 8/18 



(27r) m dt\tdt 



2m- 1 



P<1 



COs(tJiAit) ■ ■ ■ COs(^2mA 2 m^) 
VI- P IP 



xd^ 1 )(A 1 )---^ (2m) (A 2 m) 



1 8/18 



(27r) m dt\tdt 
proving (2.2). I 



2m-l r 



i2m— 1 



rfcj! • • • dw 2 m 

COs(wiAi^) • • • COs(w 2 m^2m£) 
VI- IM| 2 



(2.10) 



dcji • • • dcj 



2;)! 



3 The non-commutative case 

Let A, B be (not necessarily commuting) self-adjoint operators on the Hilbert 
space H. Then the operators A 2 , B 2 are semi-bounded. Recall the Trotter 
product formula 



-p(A 2 + B 2 ) 



lim 

m — >oo 



(3.1) 



valid for p > (or, more generally, for Re p > 0). Setting n = 2m + 1, 



A . 5 _ A 

A l — —j= i A 1 — ~ F= ) ^3 — — F= ) 

i/m i/m VTO 



A 2m~l — — F= ) — — p= 

\/m \/m 



A 



2m+l 
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we infer from (2.5) that 

= rv-^-^ [ 

j ( 47r p)-+l/2 J Q J s 

( t0J 2 B\ (tu>2m-lA 



X COS 



COS 



s 2m 
^ B 



tuiA 



(3.2) 



dtodt . 



2duJi---duJ2n 



Note that du> =- 

Vi-K+-+" 2 2 m ) 

from (3.1) and (3.2) that 



Integrating by parts as in Section 2, we see 



/ e *p cos {\J A 2 + B 2 t] dt = lim [ e *p - — 
Jo \ / m^ooj (2tt 



1 9/1S 



) m dt\t dt 



m— 1 



i2m-l 



wjH hw| m <l 



cos 



A/l-K + ---+^2m) 

The relation (3.3) suggests, at least formally, that 



dWl • • • dw 2 m 



(3.3) 



cos (\J A 2 +B 2 t] = lim — - 



1 d (Id 



/2 m- 



cos 



(27r) m dt\tdt 
(^)cos(^f)--cos(^)cos( 



/, 2 +. 



tu> 2m B 



VI - (w? + • • • + o,2 m ) 



dWl • • • dW2m 



(3.4) 

In the rest of this section we try to interpret and to justify the heuristically 
obtained formula (3.4). Note that for each m, the function 



a 



»(*) = / 

^ 2 +- 



cos 



+<<<! 



• COS 



tioi m B 



m 



duj\ ■ ■ ■ dui2m 



is a holomorphic function from the complex numbers to the space of bounded 
operators on H. Thus we "only" have to make precise the convergence in (3.4). 
Our first result is about weak convergence. 

Proposition 3.1 Let A 2 + B 2 be essentially self-adjoint on D(A 2 ) n D(B 2 ). 
For every positive integer m and h € H, define the function F m (t;A,B;h) G 
C{R l ;H) by 



F m {t-A,B;h) 

Ju>\ + 



1 d (\ d 



(27r) m dt\tdt 
( . () , ■ ~i tA 



2m-l 



i 2m tB \ 
s/m J 



(3.5) 



dwi • • • dt02mh 
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Then the functions F m (t; A, B; h) converge as vector valued temperate distribu- 
tions {i.e., in S'{R X ;H)) to cos (\M 2 + B 2 t) h. 

i t 2 
Proof. Let Rep > 0, set ^ = io~ so that Im a < 0, and let y = s. It follows 

from (3.2) and (3.3) that 

e - — e -— h=\— F m [V2s : -A 1 B-h)e- lsa ^ (3.6) 
L J V t Jo y v2s 

and, similarly, 

e -4^ = ^ jH cos (VI^T^ V2s) he—^L (3.7) 

when \J1ia is chosen to be positive for a on the negative imaginary axis. The 
left hand sides of (3.6) and (3.7) are uniformly bounded (by || h ||) for Im a > 

0, and are holomorphic for Imc > 0. Hence 



2i<r 



£ 2icrm £ 2i<rr, 
A 2 + B 2 



h is the 



Fourier transform of F m (\/2s ; A, B; h) j\f2s and \J ^e 2™ h is the Fourier 

transform of cos (^ /j4 ^^^ _ gy ^he Trotter product formula, the left hand 

sides of (3.6) converge to the left side of (3.7) for every a E i? 1 \{0}. Let 
ip(o~) £ S(R 2 ,H) be arbitrary. By the dominated convergence theorem 

[// s^_\ m , ., A dn [ I a 2 +b 2 \da 

J \[e 2i ame 2i„ m j / lj _ > y ^ e 2 «a h^^j——. 

The continuity of the inverse Fourier transform in S 1 (R 1 , H ) implies that the 
sequence F m (\/2s ; A, B; h) / \/2s converges, in the temperate distribution sense, 
to cos (V A 2 + B 2 \/2s ) \/~2s. Passing back to t — y/2s and noting that 
F m (t; A, B; h) and cos (y/A 2 + B 2 t) are even functions, we obtain the result. 
I 

One could possibly get a more precise version of Proposition 3.1 by using 
a strong version of the Trotter product formula and by estimating the conver- 
gence in an appropriate Sobolev norm. We state now a result on pointwise 
convergence, under stronger assumptions on h. 

Theorem 3.1 Let h e H be such that there exist constants C , K so that for 
every multi-index a with 2m components h £ D (A ai B a2 ■ ■ ■ A a2n - 1 B a2n j) and 

|| A ai B" 2 ■ ■ ■ A a2n ~ 1 B a2n h \\< CK^\a\\ . (3.8) 

Then F m (t; A, B; h) — > cos (V 'A 2 + B 2 t) h for \t\ < uniformly in compact 

submtervals of (--^U,-^U). 

Proof. Note that for every positive integer n, 

(A 2 + B 2 ) n = A 2£ 'B 2(1 -^ ■■■A 2£ -B 2(1 -^ . (3.9) 

£,=0,1, l<i<n 
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By (3.8) we have for h satisfying the assumptions that 

|| (A 2 + B 2 ) n h \\< 2 n CK 2n {2n)\ (3.10) 
Expanding cos (V 'A 2 + B 2 t) h in Taylor series, we see that for every N, 



cos (Va 2 +b 2 tj \^\ t2n{ * A2 + B ^ nh 



n<N 

2N+2 
~~ 2 



<2n)\ y r V / l-2t 2 X 2 

n=JV+l v / n=JV+l 

(3.11) 

or that the series converges absolutely and uniformly in compact subsets of 
|t|< 1 



s/2K ' 

Expanding similarly cos (^=f) cos ) ■■■cos {j 2 ^) h in Taylor se- 

ries, we get the estimate 

/ LUitA\ ( LU 2 tB\ f U) 2m -\tA\ f L0 2m tB\ 



X 



f UI 2 tB\ (w2 m -itA\ (uj 2m tB\ ^ 
3 — =■ • • • cos = — cos — h- } (-1) 

t 2n ^ ( Wl A) 2a ^ (uj 2 B) 2a * ■ ■ ■ K_!A) 2a --' (u; 2m B) 2a ^ h 



COS — =■ COS 



— E 



< 



m n ^ (2a)! 

|a| = n 

~ |i|2" II ^2a 1 _g2a 2 ... >1 2a 2m _ ljB 2a 2ro/l 



E — E 

mn 



m" ^ ( 2q ! 
n=AT+l „ E z^ V ' 

|a|=n 

^2n ntrtnfOnV (2m) 2 



m" (2a)! ' m™ 

n=AT+l cez^m v ' n=N+l 

|a|=n 

oo 

= C ^ (4t 2 X 2 m)™ 

n=JV+l 

so that this expansion converges absolutely and uniformly in the product of 
|| uo || < 1 with compact subsets of jtj < 2g ^— ; 2K ]y^ is positive for every fixed 
m. Hence we may perform the integration in (3.5) term by term. Evaluating 
the integral Ju^u^ , = = duo\ ■ ■ ■ duj 2m as in Section 2, we obtain the non- 

l^ll — -^/l— ||cj|| 2 

commutative version of (2.9) with at — -£= for i even, a, = -5= for i odd, i.e., 
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that 

00 (_-\\n f2n 

F m (t;A,B;h) = £ i— > — — — 

2" m" 2n — 1) • • • 1 

n=0 y ' 



(3.12) 

X 22 a\ ' 

a€Z 2m 



the series converging absolutely and uniformly in compact subsets of \t\ < 
2K 1 VTE - Noting that 2"n!(2n - 1) • ■ ■ 1 = (2n)!, we may write (3.12) in the 
form 

F m (t;A,B;h) = J2 \ , „ £ — t^ 2qi ^2q 2 ■ • ■ A 2a2m - 1 B 2am h . (3.13) 

11 = V ' cez2m 

|a|=n 

The n* 11 term (sum) in (3.11) is majorized according to the assumption (3.8) by 
l * |2n V ^CK 2n (2n)\ 



(2n)\m n a! 



a! 



so that the radius of convergence of the series in (3.11) is ^ K , and by analyticity 
of F m (t; A, B; h) (as a function of t) the series expansion (3.13) is valid for 
|t| < y^^, , uniformly in m. 

Let X , Y be non-commuting finite dimensional matrices. In this case the 
convergence in the Trotter product formula H 



< X+Y ^ = lim ( e zX/rn e zY/ m y 
n — >oo \ / 



is uniform in compact subset of the complex plane. Hence the Taylor coefficients 
of the right hand side converge to those of the left hand side, or 

(X + Y) n _ ^ x ai Y a2 ■■•I a! »- 1 7" 2 » 

n\ m—Kx ^— ' m n a\ 

aez 2lJl 



for every (fixed) n. It follows from (3.8) that for every (fixed) n and h, 

(A 2 + B 2 ) n ^ A 2a lB 2a 2 ... A 2a 2m - lB a m 

1 h = lim > ; h . (3.15) 

n\ m—>oo z — ' m n a\ 
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For any given e = and T < choose TV such that the right hand side of 

(3.11) and 

\\F m (t;A,B;h)-J2ji^— E ~^B 2a2 - ■ ■ B 2 ^h\\ 
^-^ (2n) to™ as. 



are each less than | for \t\ < T. By (3.15) there exists a constant M > such 
that for every < n < N we have 

t 2n n\A 2a ^B 2a ^ ■ ■■B 2a ^ 1 t 2n ,. 2 ,.„, „ 



11 (2n)!m" ^ m a! (2n)! v ; " 3(iV+l) 

|ot|=n 

for \t\ < T and to > M. Hence 

|| F m (t;A,B;h) - cos U/A 2 + S 2 t) ft ||< e 
for |i| < T and to > M. I 

Remark 3.3 If for every K > there exists a C > such that (3.8) holds for 
all a then F m (t; A, B; h) — ► cos (V ^4 2 + £? 2 t) /i uniformly in compact subsets of 
the complex t plane. In particular, if A and B arc bounded, then this holds for 
every h in H (uniformly in bounded sets of H). 

Remark 3.4 The Trotter product formula for any number of non-commuting 
operators is well-known |t]]. In analogy to (3.4), we may obtain the formula 



cos [JA\-\ h A 2 1 1 = lim 



q J rn^oa (2n) m l/ 2 dt\t 8t 

[ 

Ju?-\ hwj,„<l 



mq—1 



tu>iAi\ ( tw„A q \ ( tOJg+lAi 

COS I — =- • • ■ cos — cos 1 



; ! ■■■+-... 

x cos ( — jzz^ 1 - ) ; = — —dioi ■ ■ ■ d(jj mq 



"' ' ^l-^ + .-.+c^) 

(3.16) 

where to run over all positive integers if q is even, and over even m only if q 
is odd. (Note that if mq is odd then we may replace the right hand side of 
(3.16) by a term involving integration over S mq , as in (2.1).) The formal limit 
(3.16) may be interpreted, under suitable assumptions, as in Proposition 3.1 or 
in Theorem 3.2. We leave the details to the diligent reader. 

Remark 3.5 Analytic domination has been applied in the context of Trotter 
products in |lj. 



12 



4 Illustrations and examples 

We observe that in all w-integrations appearing earlier we may replace a term of 
the form cos^jT) by exp(io>.,T) (or by exp(— iu>jT)) since the rest of the inte- 
grand is even in u>j and sin(cjjT) is odd. We will use this observation constantly 
in this section in order to simplify a few calculations. 

Example 4.1 The wave equation in R n — the method of ascent. 

Recall that the operator i-p has a (unique) self-adjoint realization in L 2 (i? 1 ) 
and the distribution kernels of cos (f exp (it\-^) and exp {—it\-^) (solu- 
tions of the one-dimensional wave equation) are given by 



The simplest cases are n — 2,3. In the two dimensional case we obtain from 
(1.1) or (2.2) (m = 1) that 





(4.1) 






But 





+ U) 2 



l H(t*-{x-xy-(y-yY) 
t ^-{x-xf + {y-yf 



du>idu>2 



where H is the Hcaviside function. Hence 



cos 



( 




1 8 H (t 2 - (x -xf - (y -yf) 
2ndt ^t^-{x-xY + (y-yf 



(4.2) 
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The three dimensional case is actually simpler: 



ms I ' I '■'■■■' / - :: :i 



_Li?. 

4vr <9< 
J_d_ 



£ / <5(x — x — tuji)S(y — y — tuj 2 )S(z — z — tuj^)duj 
S 2 



t I S (x — (x + tuji)) 6 (y — (y + tuj 2 )) 5 (z — (z + tuj 3 )) du 

s 2 



The formulas (4.2) and (4.3) are classical. 

More generally, for n > 3 odd, n = 2m + 1, we have from (2.1) that 



(4.3) 



cos 



= cos t. 



\ 



1 dX i 

3 = 1 ■> 



1 d (\ d 



2(27r) m dt\tdt 



m — 1 



S 2 ' 



(e" Ml *^) (x ll x 1 ) ■ ■ ■ (e-^"^) (x n ,x n )<Lu 



(4.4) 



1 d (I d 



m — l 



t2m-l 



2(27r) m dt\tdt 

8{x\ —x\ — uj\t) ■ ■ ■ S(x n — x n — ui n t)duj 



S 2 " 



But 



27T™/ 2 

1 • 3 • • • (n - 2)u n = 1 • 3 • • • (n - 2)— ^- = 2(2^) r 



(4.5) 



and we recapture the classical formula [g, |3| . For general even dimensions n 
2m, we get from (2.2) that 



(V 



(a;, a;) 



1 9/15 



f 2m-l 



(27r) m dt\tdt 
6(xi -xi- ujxt) ■ ■ ■ 5{x 2m - x 2m - QJ 2m t) 



+W§ OT <1 



(4.6) 



dcji • ■ ■ dui 



2m 



Applying (4.5) (with n = 2m + 1) we see that the right hand side of (4.6) 
coincides with the well-known kernel of the operator mapping u(x, 0) to u{x,t) 
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(when ^r(x, 0) = 0) for even n [|. Thus, we are able to build the solution 
of the initial value problem for the n-dimensional wave equation from the (very 
elementary) one-dimensional solution. 

The Klein-Gordon operator K is given by K = A — a 2 , and we wish to 
represent the solution of the initial value problem for the Klein-Gordon equation 



d 2 u 



Ait — a 2 u 



(4.7) 



by calculating the distribution kernel of 



cos 



COS 



\U 



1 d 



Set Aj — i gf-, j = 1, • • ■ , n, A n+ \ = a. If n is even, n — 2m, we apply (2.1) to 
get 



1 9/13 



2(27r) m dt\tdt 



,2m-l 



« 2m 

X / COs(aW2m+l*) TT S( 



LUjt)dLO 



1 d [1 d 



(27r) m dt\tdt 

2m r 

><n 

3=1 



J 2m 



t 



i a/u 



(27r) m dt\tdt 



m — 1 



cos ( aiy 1 - E J= i 



2m (xj—Xj) 2 

p — 



ff(£ 2 - \x-x\ 2 ) 1 



- \ 2 



tJi E?=i (^) 



Hence 



cos (V— ifi^ (x, x) 



1 9/18 



(27r)™ dt\tdt 



(a^t 2 - |x-x| 2 ) 

y/t 2 - \X - xf 



H{t 2 - \x-x\ 



(4.8) 
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If n is odd, n = 2m — 1, we get from (2.2) that 



cos 



-Kt 



1 d (Id 



t2m— 1 



(27r) m dt\tdt 

COs(a,U>2mt) 



?+--H*3 m <l v/1 - (w? + • • • + w| m ) 



2m-l 

X IT ~ ^ _ tLOi)duJi ■ ■ ■ d,Ul2r 



i 



(27r) m dt\tdt 
Recall that for every c > 0, 



cos(atw 2 m) 



1 



\x-x\ 2 2 
~ P W 2m 



cos(atcj)(ia; 



= irJo(atc) 



where Jo is the Bessel function of the first kind of order 0. Hence 

m— 1 



COS 



-Kt 



1 



9/19 



2(2 7 r) m - 1 dt\t dt 
o (a^t 2 - \x-x\ 2 ^j H(t 2 -\x- x\ 2 ) 



(4.9) 



We may continue analytically the formulas (4.8) and (4.9) in a to get (for a 
purely imaginary) the solution of the initial value problem for the damped wave 
equation = Au + a 2 u, compare j|. 

Example 4.2 The harmonic oscillator. 

Consider the operator P = — 4-% + x 2 . We wish to represent the solution of 
the initial value problem for the second order hyperbolic equation 



a u 



dx'- 



+ x u 



(4.10) 



and for this purpose we express the distribution kernel of cos ( v? t 



d 2 

dx 2 



x 2 t I using the formula (3.4) with A 



l-f, B = x. Note 
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that cos(Bp)(x,x) — cos(xp)S(x — x). Hence (formally) 

m— 1 



COS 



dx 2 



x 2 t 



lim 



1 d (Id 



m^co (2n) m dt\tdt 



t-2/n- 



f 

JojfA Vu}\ m <\ 



5 (x\ — x 



— _ UJjt 



cos 



x 5 [ x 3 - x 2 - -^f= ) cos \^=Xi I 8{xa - x 3 ) 
Jm ) \ Jm 



<^2m-l \ ( U2 m t 
X 5 X2m-1 ~ X 2m -2 T=^ I COS | —^X 



m 



rri 



x S(x — X2m-i)duJi ■ ■ ■ du) 2m dx 1 ■ ■ ■ dx 2m -\ 



Hence 



cos 



d 2 
dx 2 



+ x 2 +t 



lim 

m— >oo 



1 9/19 
(27r) m dt [tdt 



I. 

Ju> 2 + ~- + 



COS 



U>2 



U>lt 

X + — = 
m V vm 



cos 



m— 1 r 



<^2< /_ , (Wl + W 3 )t 

r= :E + r= — 



wi + w 3 H hw 2m _i 



i — x I rfwi • • • duj 2r 



X • • • cos . . 

V v m y v v m 

(4.11) 

We may interpret (4.11) according to Proposition 3.1, as P is essentially self- 
adjoint in I^iR 1 ), (}£ and x are self-adjoint in their respective domains). 

Alternatively, set H — L 2 {B}\ e~ x dx) and let H n (x) denote the Hermite poly- 
nomial (— l) n e x -£-^e~ x ■ Then the well-known formulas H' n (x) = 2nH n _i(x), 
xH n (x) = H "+ l(a:) + nff„_i(a;) and / H 2 (x)e~ x2 dx = jH2 n n\ imply that if 
h e Sp[i7 ,#i, •••,#«], then ft e D (A ai B a2 ■ ■ ■ A" 2 ™- 1 B a2m ) for all a and 
there exist constants C and if such that (3.8) holds. Hence Theorem 3.2 is 
applicable. 

Example 4.3 Sum of squares of vector fields. 

Let X, Y be vector fields with analytic coefficients defined on an open subset 
f2 of R n such that X* = —X, Y* = —Y (the case of more than two fields may 
be treated similarly, see Remark 3.4). Let ip(t,x), i>(t,x) denote the solutions 
of the ode's 



^(t,x) = X<p(t,x), ip(0,x)=x 
^(t,x) = Ytp(t,x), ip(0,x)=x 



(4.12) 
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defined in a neighborhood of {0} x ft. We have, at least formally, the relation 



cos [\/X 2 + Y 2 t (x,x) = lim 



i ana 



i2m-l 



m-KX> (27T) m dt\t dt 



wfH hw? 



Vi-K + - + w 2 2 j 



x du>i ■ ■ ■ duj 



2m 



(4.13) 

Under further assumptions the integral in the right hand side of (4.13) will be 
defined for all t. A possible interpretation of (4.13) via Theorem 3.2 is possible 
along the lines of |J. 

Consider the case X 
hypoelliptic operator. Moreover, ip(t,x) — (a?i + t,x%), %f)(t,x) = (xi,X2 + t), 
and we get from (4.13) that 



g^, Y = X!-^. Then X 2 + Y 2 is a self-adjoint 



a 2 



a 2 



1 9/15 



cos I i I -^-a + x 2 ^—^ t I (x, x) = lim . 

IV dx\ x dx\ ) y ' m^oo (27r) m dt \ t dt, 







^ m — 1 









wfH hw? 



xi — x 



i=l ' l<i<j<m 



V / l-K + -'-+^ 2 ™) 



x ciwi • • • 



Similarly, set (in R 3 ) X 

is the Laplacian on the Heisenberg group. Here 



(4.14) 



^-+2x 2 J-, Y = — 2xi^-. Then A H = X 2 +Y 2 



<p = (xi + t, X2, x 3 + 2x2*), ip — {xi,x 2 + t, x 3 — 2xii) . 

The operators X, Y and A# are left invariant. Hence it suffices to compute 
kernels for x = 0. 
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After an elementary calculation, we get from (4.13) that 
cos(V=A^iWo) = lim -V %-(-%-Y It 2 " 1 - 1 [ 
( m t m 2t 2 ( 

^ i=l V i=l V H<i<j<m-1 l<i<j<m 

v/l-^i + ••• + <,) 

X C?CJi • • • <iu>2m • 

(4.15) 

The formulas (4.13)-(4.15) may be regarded as "Feynman-Kac" formulas. An 
explicit evaluation of (4.14) and (4.15) appears to be a challenge. 
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